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The first year data from the Wilkinson Microwave Anisotropy Probe are used to place stringent 
constraints on the topology of the Universe. We search for pairs of circles on the sky with similar 
temperature patterns along each circle. We restrict the search to back-to-back circle pairs, and to 
nearly back-to-back circle pairs, as this covers the majority of the topologies that one might hope 
to detect in a nearly flat universe. We do not find any matched circles with radius greater than 25°. 
For a wide class of models, the non-detection rules out the possibility that we live in a universe with 
topology scale smaller than 24 Gpc. 



Recently, the Wilkinson Microwave Anisotropy Probe 
(WMAP) has produced a high resolution, low noise map 
of the temperature fluctuations in the Cosmic Microwave 
Background (CMB) radiation pj. Our goal is to place 
constraints on the topology of the Universe by search- 
ing for matched circles in this map. Here we report the 
results of a directed search for the most probable topolo- 
gies. The results of the full search, and additional details 
about our methodology, will be the subject of a subse- 
quent publication. 

The question we seek to address can be plainly stated: 
Is the Universe finite or infinite? Or, more precisely, 
what is the shape of space? Our technique for studying 
the shape, or topology, of the Universe is based on the 
simple observation that if we live in a Universe that is 
finite, light from a distant object will be able to reach us 
along more than one path. The one caveat is that the 
light must have sufficient time to reach us from multiple 
directions, or put another way, that the Universe is suf- 
ficiently small. The idea that space might be curled up 
in some complicated fashion has a long history. In 1900 
Schwarzschild considered the possibility that space 
may have non-trivial topology, and used the multiple im- 
age idea to place lower bounds on the size of the Universe. 
Recent progress is summarized in Ref. 0- 

The results from the WMAP experiment [l| have deep- 
ened interest in the possibility of a finite universe. Several 
reported large scale anomalies are all potential signatures 
of a finite universe: the lack of large angle fluctuations Q , 



reported non-Gaussian features in the maps 0.0- and 
features in the power spectrum Q. 

The local geometry of space constrains, but does not 
dictate the topology of space. A host of astronomical 
observations support the idea that space is locally homo- 
geneous and isotropic, so we may restrict our attention 
to the three dimensional spaces of constant curvature: 
Euclidean space E 3 ; Hyperbolic space H 3 ; and Spheri- 
cal space S 3 . A useful way to view non-trivial topologies 
with these local geometries is to imagine the space be- 
ing tiled by identical copies of a fundamental cell. For 
example, Euclidean space can be tiled by cubes, result- 
ing in a three-torus topology. Assuming that light has 
sufficient time to cross the fundamental cell, an observer 
would see multiple copies of a single astronomical ob- 
ject. To have the best chance of seeing "around the Uni- 
verse" we should look for multiple images of the most 
distant reaches of the Universe. The last scattering sur- 
face, or decoupling surface, marks the edge of the visible 
Universe. Thus, looking for multiple imaging of the last 
scattering surface is a powerful way to look for non-trivial 
topology. 

But how can we tell that the essentially random pat- 
tern of hot and cold spots on the last scattering surface 
have been multiply imaged? First of all, the microwave 
photons seen by an observer have all been traveling at 
the same speed, for the same amount of time, so the 
the surface of last scatter is 2-sphere centered on the ob- 
server. Each copy of the observer will come with a copy 
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of the surface of last scatter, and if the copies are sepa- 
rated by a distance less than the diameter of the surface 
of last scatter, then the copies of the surface of last scat- 
ter will intersect. Since the intersection of two 2-sphcrcs 
defines a circle, the surfaces of last scatter will intersect 
along circles. These circles are visible by both copies of 
the observer, but from opposite sides. Of course the two 
copies are really one observer, so if space is sufficiently 
small, the cosmic microwave background radiation from 
the surface of last scatter will have patterns of hot and 
cold spots that match around circles @- The key as- 
sumption in this analysis is that the CMB fluctuations 
come primarily from the surface-of-last scatter and are 
due to density and potential terms at the surface of last 
scatter. 

Implementing our matched circle test is straightfor- 
ward but computationally intensive. The general search 
must explore a six dimensional parameter space. The pa- 
rameters are the location of the first circle center, (6*i ,<j>i), 
the location of the second circle center, (02, <h), the an- 
gular radius of the circle a, and the relative phase of the 
two circles j3. We use the Healpix scheme to define our 
search grid on the sky. A resolution r Healpix grid di- 
vides the sky into N = 12N 2 ide equal area pixels, where 
N s ide — 2 r . Similar angular resolutions were used to step 
through a and (3. Thus, the total number of circles be- 
ing compared scales as ~ TV 3 , and each comparison takes 
~ N 1 / 2 operations. The simplest implementation of the 
search at the resolution of the WMAP data (r = 9) would 
take greater than 10 20 operations. This is not computa- 
tionally feasible. However, with the algorithms outlined 
below, we can carry out a complete search for most likely 
topologies using the WMAP data. 

The WMAP data suggest that the Universe is very 
nearly spatially flat, with a density parameter J7o = 
1.02 ± 0.020. Our universe is either Euclidean, or its 
radius of curvature is large compared to radius of the 
surface of last scatter. For topology to be observable 
using our matched circle technique we require that the 
distance to our nearest copy is less than the diameter 
of the surface of last scatter, which in turn implies that, 
near our location and in at least one direction, the funda- 
mental cell is small compared to the radius of curvature. 
Given the observational constraint on the curvature ra- 
dius, it is highly unlikely that there are any hyperbolic 
topologies small enough to be detectable 0, and there 
are strong constraints on the types of spherical topologies 
that might be detected 01- Naturally, the near flatness 
of the Universe does not place any restrictions on the 
observability of the Euclidean topologies. Remarkably, 
the largest matching circles in most of the topologies we 
might hope to detect will be back-to-back on the sky 
or nearly so. This immediately reduces the search space 
from six to four dimensions. This result is exact and easy 
to prove for nine of the ten Euclidean topologies. The one 
exception is the Hantzsche-Wendt space, which has its 
largest circles with centers separated by 90°. The result 
is less obvious in spherical space |l0j but is nonetheless 



exact for a large class of such spaces (single action man- 
ifolds). All others (double or linked action manifolds) 
predict slightly less than 180° separations between the 
circle centers for generic locations of the observer, but 
larger offsets for special locations. This generically slight 
offset motivates the final search described in this letter. 
The possibility of a large offset requires the full search 
which will be described in a future publication. 

How do we compare the circles? Around each pixel 
i, we draw a circle of radius a and linearly interpolate 
values at n = 2 r+1 points along the circle. We then 
Fourier transform each circle: Ti(<p) — Ti m exp imtp 
and compare circle pairs, with equal weight for all angular 
scales: 
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The i, j refer to the location of the circle centers. With 
this definition Sij = 1 for a perfect match. For two ran- 
dom circles < Sij >= 0. We estimate < Sfj > below. We 
can speed the calculation of the S statistic by rewriting 
it as 
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By performing an inverse fast Fourier transform of the 
s m , we get Sij (a, (3) at a cost of ~ N 1 / 2 logiV operations. 
This reduces the cost of the back-to-back search from 
~ TV 5 / 2 to ~ iV 2 logiV. We further speed the calculation 
by using a hierarchical approach: we first search at r = 7 
and identify the 5000 best circles at a fixed value of a. 
We then refine the search around these best pixel centers 
and radii and complete the search at r — 9. We store the 
best 5000 matches at each a. S ma x(o) is the best match 
value found at each radius. 

To test the search algorithm, we generate a simulated 
CMB sky for a finite universe. When only the Sachs- 
Wolfe effect was included, the algorithm found nearly 
perfect circles: S rnax — 0.99 with pixelization effects 
accounting for the residual errors. Real world effects, 
however, significantly degrade the quality of the circles. 
The dominant source of noise are velocities at the surface 
of last scatter which dilutes the quality of the expected 
match. Following the approach outlined in the appendix 
of ^2 ' we generate CMB skies with initial Gaussian ran- 
dom amplitudes for each k mode and then compute the 
transfer function for each value of k. This approach im- 
plicitly includes all of the key physics: finite surface of 
last scatter, reionization, the ISW effect and the Dopplcr 
term. Since topology has little effect on the power spec- 
trum for multipoles I > 20, we use the best fit param- 
eters based on the analysis of the WMAP temperature 
and polarization data Q| . The ratio of the amplitude of 
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FIG. 1: The maximum value of the circle statistic as a func- 
tion of radius, a, for a simulated finite universe model. The 
peaks in the plots correspond to positions of matched circles. 
The cyan line corresponds to the detection threshold discussed 
below, a is measured in degrees in all figures. S m ax{a) is the 
best match value found at each radius. 
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FIG. 2: The maximum value of the circle statistic as a func- 
tion of radius for back-to-back circle pairs on scrambled CMB 
skies. The results are shown for seven different realizations. 
The red line is a fit for the false positive rate. The cyan line 
shows the false detection threshold set so that fewer than 1 
in 100 simulations would yield a false event. 



the potential and density terms (which contribute to the 
signal for the circle statistic) to the velocity and ISW 
terms (which contribute to the noise for the circle statis- 
tic) depends weakly on the basic cosmological parame- 
ters: within the range consistent with the WMAP data, 
the variation is small. For the simulated sky, we assume 
that the primordial potential fluctuations had the usual 
1/k 3 power spectrum. The simulation included detec- 
tor noise and the effects of the WMAP beams ^3 ■ 
The simulation did not include gravitational lensing of 
the CMB as the len sing deflection angle in the standard 
cosmology is small 15j, ~ 6', compared to the smooth- 
ing scale used in the analysis, 40'. When the search was 
performed on these realistic skies, the quality of circle 
matches were degraded. Figure 1 shows the results of the 
search on the simulated sky: the peaks in the plots cor- 
respond to radii at which a matched circle was detected. 
The largest circles had S max — 0.75, and the best match 
value decreased for smaller circles. The code was able to 
find all 107 predicted circle pairs with radii greater than 
25°: the poorest detection had an S max value of 0.45. 
While we only simulated the three torus, the degrada- 
tion of the circle match as a function of radius will be 
nearly identical for other topologies as the contribution 
of the Doppler term (and detector noise) is determined 
by local physics at angular scales well below the circle 
radius. 

For the circle search, we generated a low foreground 
CMB map from the WMAP data. Outside the WMAP 
Kp2 cut, we used a noise- weighted combination of the 
WMAP Q, V and W band maps. Using the template 
approach of these maps were corrected for dust, 
foreground and synchrotron emission. This map was 



smoothed to the resolution of the Q band map. Inside 
the Kp2 cut, we used the internal linear combination 
map 16]. 

What is the expected level of the false positive sig- 
nal? If we approximate the CMB sky as a Gaussian ran- 
dom field with coherence angle 9 C , then there are roughly 
Ndrc = 27r/(9 2 independent circles on the sky of radius 
a. Along each circle, there are 27rsin(a)/0 c independent 
patches and also 27rsin(a)/# c independent orientations. 
Thus, the back-to-back search considered N searc h — 
47T 2 sin(a) / 6% possible circle pairs. If we treat the patches 
as independent, then < S 2 >V 2 = l/^27r sm(a)/9 c . The 
number of circles expected above a given threshold is, 
N(S > S ) = N search /2erfc(S Q /V2 < S 2 >). Thus, in- 
verting the complementary error function yields the max- 
imum value for any false positive circle at each radius, 
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For both the simulations and the data, C ~ 0.7°. Substi- 
tuting the expression for N searc h into equation J2J yields: 
^max ~ 0.24-^/1 + ln(sina)/18/\/sina, This estimate is 
sensitive to 9 C and excludes correlations between circles 
We use "scrambled" versions of the WMAP cleaned 
sky map to obtain a more accurate estimate of the false 
positive rate. We generate "scrambled" versions of the 
true sky by taking the spherical harmonic transform of 
the map and then randomly exchanging a\ m values at 
fixed I. This scrambling generates new maps with the 
same two point function but different phase correlations. 
The S max value for the best fit back-to-back circle found 
at each radius is plotted as a function of a in figure 2. 
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FIG. 3: The maximum value of the circle statistic as a func- 
tion of radius for the WMAP data for back-to-back circles. 
The solid line in the figure is for a orientable topology. The 
red line is for non-oriented topologies. The green line is the 
expected false detection level. The cyan line is the detection 
threshold. The spike at 90° is due to a trivial match between 
a circle of radius 90° centered around a point and a copy of 
the same circle centered around its antipodal point. 




FIG. 4: The maximum value of the circle statistic as a func- 
tion of radius for the WMAP data for circles whose centers are 
separated by greater than 170°. The solid line in the figure is 
for a orientable topology. The blue line is the expected value 
of S m ax f° r a Gaussian sky (based on the analytic estimate 
and extrapolation from the scrambled sky simulations). The 
cyan line is the detection threshold. 



The false positive signal is well fit by a slightly modified 
form of the analytical estimates. Based on the analytical 
estimate, we set a detection threshold so that fewer than 
1 in 100 random skies should generate a false match. This 
threshold is shown in figures (1), (2) and (3). 

Figure 3 shows the result of the back-to-back search 
on the WMAP data. The search checked circle pairs 
that were both matched in phase (points along both cir- 
cle labeled in a clockwise direction) and nipped in phase. 
These two cases correspond to searches for non-orientable 
and orientable topologies. We did not detect any statis- 
tically significant circle matches in cither search. The 
simulations show that the minimum signal expected for 
a circle of radius greater than 25° is S max = 0.45, which 
is above the false detection threshold. Since we did not 
detect this signal, we can exclude any universe that pre- 
dicts circles larger than this critical size. In a broad class 
of topologies, this constrains the minimum translation 
distance, d: 

d = 2i? c atan [tan(i?/ ss /i? c ) cos(a mi „)] (5) 

where R c is the radius of curvature, and Ri ss is the dis- 
tance to the last scattering surface. Note that one recov- 
ers the Euclidean formula for R c — > oo and the hyperbolic 
formula for R c —> iR c . For a flat or nearly flat universe, 
models with d < 2Ri ss cos(a) ~ 24 Gpc are excluded. 
Note that the search excludes the Poincare Dodecohe- 
dron suggested in Q as this model predicts back-to-back 
circles of radius 35°. 

By considering a broader set of circles, we can extend 
the search to other possible topologies. Figure 4 shows 



the result of a search for nearly back-to-back circles: cir- 
cles whose centers are separated by more than 170°. Even 
this restricted search requires searching far more circles 
than the back-to-back search: the search takes 4 proces- 
sor months on a SGI Origin computer. Because the num- 
ber of independent circles in this search is roughly 200 
times larger than in the back-to-back search, the level 
of the false positive signal is higher. The increase in 
N searc h raises the false detection threshold from 5.41 to 
6.29 < S 2 > x ' 2 : S^ ax = 0.2 Vl + ln(sina)/6/\/shia. 
Because of this higher signal, the minimum detectable 
angle is 28° for this topology. The signal at a ~ 85 — 90° 
is due to nearly overlapping and osculating circles. This 
feature is also seen in scrambled skies. 

Has this search ruled out the possibility that we live 
in a finite universe? No, it has only ruled out a broad 
class of finite universe models smaller than a character- 
istic size. By extending the search to all possible orien- 
tations, we will be able to exclude the possibility that we 
live in a universe smaller than 24 Gpc in diameter. More 
directed searches could extend the result for specific man- 
ifolds somewhat beyond 24 Gpc. With lower noise and 
higher resolution CMB maps (from WMAP 's extended 
mission and from Planck), we will be able to search for 
smaller circles and extend the limit to ~ 28 Gpc. If the 
universe is larger than this, the circle statistic will not be 
able to constrain it shape. 
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